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Abstract. In this paper, we study the behaviour of the output of pure entangled states after 
being transformed by a product of conjugate random unitary channels. This study is motivated 
by the counterexamples by Hastings [25] and Hayden- Winter [27] to the additivity problems. 
In particular, we study in depth the difference of behaviour between random unitary channels 
and generic random channels. In the case where the number of unitary operators is fixed, we 
compute the limiting eigenvalues of the output states. In the case where the number of unitary 
operators grows linearly with the dimension of the input space, we show that the eigenvalue 
distribution converges to a limiting shape that we characterize with free probability tools. In 
order to perform the required computations, we need a systematic way of dealing with moment 
problems for random matrices whose blocks are i.i.d. Haar distributed unitary operators. This 
is achieved by extending the graphical Weingarten calculus introduced in [T3] . 

< 

1. Introduction 

In Quantum Information Theory, random unitary channels are completely positive, trace 
preserving and unit preserving maps that can be written as 
.*_> ■ /,. 

$(p) = J2 Wi U iP U*, 

i=l 

where Ui are unitary operators acting on C n and Wi are positive weights that sum up to one. 
This class of quantum channels is very important, not only because the action of such a 
Q"\ | channel can be understood as randomly applying one of the unitary transformations U with 

respective probabilities Wi, but also because it has highly non-classical properties. For example, 
random unitary channels have been used to disprove the additivity of minimum output entropy 

ESI. 
oo ■ — 

In this paper, we are interested in a natural setting in which we take a convex combination 
£N| . of k random unitary evolutions; in other words, we choose the unitary operators U at random 

and independently from the unitary group, following the Haar distribution. The behaviour of 
this kind of quantum channels has been extensively studied in the literature [TJ \*ZU[ |2"5"] . 

We are principally interested in the study of the moments of typical outputs for products of 
conjugated random unitary channels. One of the main results of our paper is that the typical 
outputs are deterministic when one takes a product of such a channel and its complex conjugate 
and applies it to entangled input states, in the spirit of [141 1161 fT^j ■ More precisely, our main 
results can be stated informally as follows (for precise statements, see Theorems 15.21 and 16. 7|) : 

Theorem 1.1. Consider the output state Z n = [3> <%> ^](ip n tfj^) , given as the image of a "well 
behaved" pure state ijj n under the product of conjugate random unitary channels. 

If k, the number of unitary operators Ui, is fixed and n — > oo, then the set of the k 2 non-zero 
eigenvalues of Z n is 

{wiWj : i,j = 1, ...,k,i / j}U{«j : i = l,...,k}, 

where the numbers Sj depend on w and some parameter m quantifying the overlap between the 
input state ip n and the Bell state tp n . 
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// both dimensions k and n grow to infinity, but at a constant ratio k/n — )• c, then the largest 
eigenvalue of Z n behaves as Const. /k and the rest of the spectrum, when properly rescaled, 
converges towards a compound free Poisson distribution that can be characterized in terms of 
the behaviour of the ratio c and the weights Wi . 

We strengthen previous results on the eigenvalues of the output states Z n , obtained by linear 
algebra techniques, by computing exactly the asymptotic spectrum of the output states. Such 
precise results, along with results for single channels |13j . may improve existing violations of the 
additivity of the minimum output entropy or similar quantities. For the purpose of studying 
these random channels, we extend the graphical calculus introduced in [T3] to the case of multiple 
independent, identically distributed ( i.i.d. ) random Haar unitary matrices, see Theorem 13. Ill 

Our paper is organised as follows. In Section [2] we recall facts about quantum channels and 
give the definition of random unitary channels. In Section [3l we extend Weingarten calculus to 
include the case of several independent Haar unitary matrices. Section H] provides a comparison 
of Haar random isometries and block random isometries. Even though it is intuitively clear 
that the two are quite different from each other, we give explicit instances of the differences via 
a few explicit moment computations. These moment computations serve also as a warm-up for 
Sections [5] and O where the behaviour of output of entangled pure states is considered under 
different scalings; these sections contain the main results of the paper. In SectionOwe treat the 
case of a fixed number of unitary operators, while in Section [6] we consider the scaling k/n — > c. 
We analyze thoroughly the entropies of the limiting objects, finding the parameter values which 
yield the output states with the least entropy. Finally, we make a few concluding remarks in 
Section 

2. Random unitary channels 

2.1. Our model. Random unitary channels are natural models among quantum channels and 
have played a key role in the research around violations of different quantities related to classical 
capacities of quantum channels. 

Random unitary channels are unit-preserving quantum channels that can be written as 



where Ui 6 lA{n) are independent random unitary operators with uniform (or Haar) distribution 
and Wi are fixed positive weights that sum up to one. Hastings used a similar model to prove 
additivity violation of minimal output entropy [25| . In our model, a unitary operator U{ is 
applied to an input with some fixed probability Wi, while in Hastings' model those probabilities 
were chosen randomly. 

2.2. Stinespring picture and complementary channel. The above random channel 3> : 
M n (C) — > M n (C) is described in the Stinespring picture as 



k 



(1) 




i=l 



(2) 
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where 
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is an isometry, 
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is a stack of unitary operators and Tr^ is the partial trace operator Tr<S>id acting on M^(C) <8> 
M n (C). In matrix representation, 

/ wiUtpU* ... y/wiWkUipU^ 

(4) VpV* = : • • . : 

\y/w k wiUkpU* . . . WkUkpUl 

As one sees, tracing out the first space (the environment space C k ) gives the sum of diagonal 
blocks, which is the quantum channel $. On the other hand, if we switch the roles of output 
and environment spaces, i.e., if we trace out the second space (the output space C n ), we get the 
complementary channel $ c : M n (C) ->• M k (C) [29~1 13T] : 

(5) ($ c (p)) i;j = ^Wiw]Tr[U iP U*] 

A quantum channel and its complementary channel share the same output eigenvalues, up to 
zeroes, for any pure input via Schmidt decomposition, see [291 EI] for details. 

Recall that the Shannon entropy functional, defined for probability vectors p = (pi) 

H (p) = ~ ^Pdogpi, 

i 

extends to unit trace, positive matrices (or quantum states) p via functional calculus 

H(p) = -Tr(plogp), 
For a quantum channel <E>, one can define its minimal output entropy 

# min ($) = min H(*(p)). 

p state 

By using convexity properties, one can show that the minimum above is reached on pure 
states (i.e. rank one projectors) so that the value of the minimum output entropy does not 
change when by replacing $ by <& c . Our interest in complementary channels is motivated by 
the fact that, often, the size of the environment is smaller than the output size, so output states 
are easier to study for Q c . 

In the seminal paper of Hastings [25], violation of additivity was found when the dimensions 
of input and output spaces are much larger than that of the environment space of the channel. 
In this paper, we shall study two asymptotic regimes, k fixed and n — > oo and then k, n — )■ oo, 
with a linear growth k/n — > c. 

3. WeINGARTEN CALCULUS FOR SEVERAL INDEPENDENT UNITARY MATRICES 

The method of graphical calculus was introduced in [14] and later used in [161 EH OS H2] 
to study random matrix models in which Haar distributed unitary operators played a major 
role. In particular, the limiting eigenvalue distribution of [<J> ® ^](ip n ip^) , where $ is random 
quantum chanel and ^ n is a generalised Bell state, was calculated in |12j . The aim of this paper 
is to perform similar computations for random unitary channels instead of random channels. 

In the previous work in which random quantum channels were analysed using the graphical 
calculus method |14t LT6l \T7\ [12] , the isometry defining the channel in the Stinespring picture 
was a truncation of a Haar-distributed random unitary matrix. Hence, only one random unitary 
operator was needed to perform the computations. In this paper, we need to develop a new 
technique of graphical calculus for k i.i.d. Haar-distributed unitary matrices. 

We start with the usual, 1-matrix, graphical Weingarten calculus and then generalise it to 
cover the case of block-isometries built up from independent Haar-distributed unitary matrices. 

3.1. Weingarten formula for a single matrix. Let us start by recalling the definition of a 
combinatorial object of interest, the unitary Weingarten function. 

Definition 3.1. The unitary Weingarten function Wg(n,<r) is a function of a dimension pa- 
rameter n and of a permutation a in the symmetric group S p . It is the inverse of the function 
a i — y ri& u under the convolution for the symmetric group (j^a denotes the number of cycles of 
the permutation a). 
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Notice that the function a i— > ri& a is invertible when n is large, as it behaves like n p 5 e 
as n — y oo. Actually, if n < p the function is not invertible any more, but we can keep 
this definition by taking the pseudo inverse and the theorems below will still hold true (we 
refer to [19] for historical references and further details). We shall use the shorthand notation 
Wg(ff) = Wg(n, a) when the dimension parameter n is clear from context. 

The function Wg is used to compute integrals with respect to the Haar measure on the 
unitary group (we shall denote by U(n) the unitary group acting on an n-dimensional Hilbert 
space). The first theorem is as follows: 

Theorem 3.2. Let n be a positive integer and i = (i±, . . . ,i p ), %' = (i'i, ■ ■ ■ ,i' p ), j = (ii, • • • ,j P ), 
f = . . . , j' p ) hep-tuples of positive integers from [n] = {1, 2, . . . , n}. Then 



(6) / U hjl ■ ■ ■ Ui j U e f ■ ■ ■ U V jl dU = 

JU{n) 

E S hU) ■ ■ ■ 6 ^j^r W ■ ■ ■ Wg(n ' T<T ~ l) - 

If P ¥" P* then 

(7) / u hjl ...u ipj jf^...wy f du = o. 

JU(nj 1 ' 

Since we shall perform integration over large unitary groups, we are interested in the values 
of the Weingarten function in the limit n — > oo. The following result encloses all the data we 
need for our computations about the asymptotics of the Wg function; see [11] for a proof. 

Theorem 3.3. For a permutation a G S v , we write c E a to denote the fact that c is a cycle of 
a. Then 

(8) Wg(n,a) = [JWg(n,c)(l + 0(n- 2 )) 

C£(7 

and 

(9) Wg(n,(l,...,d)) = (-l) d - 1 c (i _ 1 [] (n-j)- 1 

where Cj = ., is the i-th Catalan number. 

As a shorthand for the quantities in Theorem 13. 3\ we introduce the function Mob on the 
symmetric group. Mob is invariant under conjugation and multiplicative over the cycles; further, 
it satisfies for any permutation a G S p : 

(10) Wg(n, a) = n-fr + l CT l)(Mob((T) + 0(n~ 2 )) 

where |tr| = p — jfo is the length of a, i.e. the minimal number of transpositions that multiply 
to cr; we shall also use the notation | • | for the cardinality of sets. We refer to [19] for details 
about the function Mob. We shall make use of the following explicit formulas. 

Lemma 3.4. The Mobius function is multiplicative with respect to the cycle structure of per- 
mutations 

Mob( C r) = n(-l)l 6 l- 1 Cat| fehl , 

beer 

where Cat n is the n-th Catalan number. In particular, if a is a product of disjoint transpositions, 
then 

(11) Mob(a) = (-1) H . 

We finish this section by a well known lemma which we will use several times in the proofs 
of the main theorems of this paper. This result is contained in [34] . 

Lemma 3.5. The function d(a,r) = |cj _1 t| is an integer valued distance on S p . Besides, it has 
the following properties: 
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• the diameter of S p is p — 1; 

• d(-,-) is left and right translation invariant; 

• for three permutations a%, 02, r G S p , the quantity d(r, a\) + d(r, 02) has the same parity 
as d(ai,a 2 ); 

• the set of geodesic points between the identity permutation id and some permutation 
g G S p is in bijection with the set of non-crossing partitions smaller than ir, where the 
partition it encodes the cycle structure of a. Moreover, the preceding bijection preserves 
the lattice structure. 

3.2. Weingarten formulas for several independent unitary matrices. In this paper we 
want to treat random matrix models in which several i.i.d. Haar random unitary matrices 
appear. We shall use the independence property and the Weingarten formula together and 
unify these properties into an unique statement, as in Theorems 13.71 and 13.81 Before we can do 
this, we introduce the following notation (we put [n] = {1, 2, . . . , n}): 

Definition 3.6. Consider two integer functions I : \p] — > [k] and I' : [p'] — > [k]. Whenever 
p = p' , we denote by Sp"* 1 ' the set of permutations of p objects which map the level sets of I to 
the level sets of I' : 

Sf 1 ' = {o£S p \o-{r\s)) = o-{l'-\s))}. 
Note that this set is empty iff 

\{t:l t = s}\^\{t:l' t = s}\ 
for some s G [k]. If p ^ p' , we put St~* 1 ' = 0. 

Note that permutations a G S p ~^ 1 ' admit a decomposition along the level sets of /' 
(12) a = a\ X • • • x crfc, 

where 

a s : r\s) -> l'-\s). 

The main result of this section is the following generalization of the Weingarten formula for 
moments in several independent Haar unitary matrices. 

Theorem 3.7 (generalized Weingarten formula). Let n and k be positive integers and i = 
(h, . . . , i p ), i' = . . . , i' pl ), j = (ji, . . . ,j p ), j' = (J[, . . . ,j' pl ) be tuples of positive integers from 
[n] and I = (Ji, . . . , l p ), I' = (l[, . . . , l' pl ) be tuples of positive integers from [k]. Then 



(13) f U? 1 } • • • Ul 1 *} UP, • • • U { H dZ7« • • • dU^ 

( 14 ) = E s *k (1) ■ ■ ■ ■ ■ ■ Wg ^"( n ' ®> 

where the modified Weingarten function Wg ! ^' is defined via the product formula 

k 

Wg M '(n,a,i3) = nWgKa; 1 /?,). 

s=l 

Proof. We start by factoring the integral using the independence of the random unitary matrices 
U(s). 



(15) / uf 1 } ■ ■ ■ uf p ] ufl ■ ■ ■ uf p l duw . . . duw 



U{n) 



= f [factors with l t = l,l' t = 1] dU {1) x • • • x f [factors with l t = k,l' t = k] dU {k) 

JU(n) ' ' JU(n) 

The above product vanishes whenever |/ _1 (s)| ^ \l' (s)\ for some s G [k]. The value of each 
factor is computed using the usual Weingarten formula in Theorem l3.2l and the result follows. □ 



6 BENOIT COLLINS, MOTOHISA FUKUDA, AND ION NECHITA 

One can reformulate the above result in a more practical way, eliminating the restriction on 
the sum indices a, f3, as follows. 

Theorem 3.8 (generalized Weingarten formula, second take). Let n and k be positive integers 
and i = (i x , . . . , i p ), i' = i' pl ), j = (ji, . . . ,j p ), f = (J[,... ,j' p ,) be tuples of positive 

integers from [n] and I = (li, . . . ,l p ), I' = (l[, . . . ,l' pl ) be tuples of positive integers from [k]. 
Then 



f i fi\ / rr( z i) . . . rrft>) jjQ'i) . . . jjV'p) (jjj. . . . AT j, 

^ J u{n)h u h,h u iv,j P u i' 1 ,3[ u i' P ,r P aUl dUk 



(*) (**) 
Here, 



Wg^a" 1 /?) = jjwg(n,a7V-) when S p ~* 1 ' , 

s=l 

and it can have any other value when one of a or (3 is not in S p ~^ 1 ' . 

Proof. This is a direct consequence of Theorem 13.71 Indeed, because of the factors (**) in the 
above formula, only a, f3 € S p survive, and the value of the Weingarten weight is the same 
as in Theorem 13.71 □ 

For the concrete applications that follow, it is useful to have a simple equivalent for the 
quantity Wg(n,a _1 /3). 

Proposition 3.9. Assuming that p is fixed, for every I, I' and a, (3 6 St^ 1 , we have 

(17) Wg^'(n,a,/3) =Wg(n,a -1 y 5) = Wg(n, a' 1 + 0(n~ 2 )) 

Proof. We shall use the fact that one can approximate the usual Weingarten function by a 
polynomial times the Mobius function, which are multiplicative. Start from the left hand side 
of the above equality, and use the fact that both a and (3 decompose as products of smaller 
permutations acting on the level sets of I and I', see (|12p : 



Wg M '(n,Q^) = JIWg(n,a 8 ,ft 



8=1 

k 

= 11 n -\rH')\-\* l P.\ Mob(a7 1 /3,)(l + 0(n- 2 )) 

8=1 

= n - p - |o ~ 1/31 Mob(a~ 1 /3)(l + 0{n- 2 )) 
= Wg(n,a- 1 /?)(l + Ort. 

Note that the product over k above has at most p non trivial factors, because at most p of 
the permutations a s ,f3 s are non trivial. 

□ 



3.3. Graphical Weingarten calculus for several independent unitary matrices. In this 
section we extend the graphical formalism introduced in [H] to encompass integrals over several 
independent unitary Haar-distributed matrices. We first recall briefly the single-matrix case. 
For more details on this method, we refer the reader to the paper [TJ] and to other work which 
make use of this technique [TBI IT71 [TBI fT2] . 

In the graphical calculus, matrices (or, more generally, tensors) are represented by boxes. 
To each box, one attaches symbols of different shapes, corresponding to vector spaces. The 
symbols can be empty (white) or filled (black), corresponding to spaces or their duals. 
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Wires connect these symbols, each wire corresponding to a tensor contraction V x V* — > C. 
A diagram is a collection of such boxes and wires and corresponds to an element in a tensor 
product space. 

The main advantage of such a representation is that it provides an efficient way of computing 
expectation values of such tensors when some (or all) of the boxes are random unitary matrices. 
In p3], the authors proposed an efficient way to apply the Weingarten formula ((6]) to a diagram, 
which has a nice graphical interpretation. 

The delta functions in each summand in the RHS of © describes how one should connect 
boxes. Each pair of permutations (a, f3) in ([6]) connects the labels of U and U boxes and then, 
after erasing those boxes and keeping the wires, one obtains a new diagram. The permutation 
a is used to connect the white (or empty) labels of the boxes in the following way: the white 
decorations of the i-th U box are connected to the corresponding white decorations of the a(i)-th 
U box. Permutation f3 is used in a similar manner to connect the black (or filled) decorations. 

This process for a fixed permutation pair (a, 0) is called a removal and the whole process 
which sums all the new graphs over the all permutations is called the graph expansion. Impor- 
tantly, the graphical calculus works linearly and separated components are related by tensor 
products, as is assumed implicitly above. The above procedure is summarized in the following 
important result. 

Theorem 3.10. If T> is a diagram containing boxes U,U corresponding to a Haar- distributed 
random unitary matrix U 6 U{n), the expectation value ofT> with respect to U can be decomposed 
as a sum of weighted diagrams T> a ^ obtained by connecting the white labels of the U boxed along 
the permutation a and the black labels along (3. The weights of the diagrams are given by the 
Weingarten functions. 

E [/ (P) = ^D^Wg( n ,ar 1 ). 

We now extend the graphical calculus by adding a new box corresponding to a column-block 
matrix made of independent unitary matrices and generalizing the removal procedure in order 
to allow to take expectations of diagrams containing such boxes. 

Consider i.i.d. Haar unitary matrices U^ l \ U^ 2 \ . . . , E U(n) and stack them up into a 
block-column V G Mfc nxn (C) 

k 

i=i 

where {e^} is an orthonormal basis of C k used to index the unitary blocks. Note that, up to a 
constant, V is an isometry from C n — > C k ®C n , i.e. V*V = k\ n . Graphically, one can represent 



V as in Figure 1(a) With this correspondence, one can read the Weingarten formulas (|13p and 



(|16p in terms of the matrix V by using the identification ujj = Moreover, the diagram 

of the (true) isometry V used to defined random unitary channels can be easily obtained from 
V and the diagonal matrix 

(18) W = diag^x,. . . ,w k ), 



see Figure 1(b) 



Consider now a diagram T> containing a box V corresponding to a stack of i.i.d. Haar unitary 
operators. We are interested in computing the expectation value of T> with respect to V. 
In the spirit of Theorem 13.10} we are going to express this expectation value as a sum over 
diagrams obtained by a new removal procedure, weighted by Weingarten coefficients. One of 
the main differences between the Weingarten formula for a single Haar unitary matrix ([6|) and 
the Weingarten formula for several such Haar random matrices (|16p is the fact that the indices 
I and V in the latter equation are connected by both permutations a and (3. Hence, we need 
to redefine the removal procedure when computing expectation values with respect to V boxes. 
Since the 1,1' indices are connected using both permutations, we shall add a duplication symbol 



inside the box of V of V"*, next to the label corresponding to C , see Figure 1(c) During the 
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V 



V 



w 



V 



V 



(a) 



(b) 



(c) 



Figure 1. Boxes for block isometries. Round labels correspond to C n and the 
square label corresponds to C k . On the left, the usual graphical representation 
of a stack of unitary operators. In the middle, the box V for a weighted block- 
isometry. On the right, the box for a stack of i.i.d. random unitary matrices. The 
square label has a duplication symbol associated to it on the inside of the box, 
suggesting it should be connected using both permutations when performing the 
graph expansion. 



removal procedure, labels associated to the duplication symbols should be connected with both 
permutations a and /3, as in Figure [2j 



V 



- r 



a(i) 



b 



Figure 2. Removal procedure for boxes corresponding to stacks of i.i.d. Haar 
unitary operators. Labels with duplication symbols must be connected using 
both permutations. 



After the removal procedure, the duplication symbols remaining must be interpreted in terms 
of the duplication operator T : C k — > C k (g> C k and its adjoint T* (see Figure [3]): 

k 



1=1 

A- 



T* = £( e 




Figure 3. Duplication operators T and T* . 



Theorem 3.11. Let D be a diagram containing boxes V, V corresponding to a random matrices 
having a block-Haar distribution. The expectation value ofV with respect to V can be decomposed 
as a sum of weighted diagrams T> a a obtained by connecting the white labels of the V boxed along 
the permutation a, the black labels along f3, and the labels with a duplication symbol along both 
a and /3. The weights of the diagrams are given by modified Weingarten functions. 

Ev(P) = ^P Qi/3 Wg(n,a,/3). 
Several applications of this result shall be discussed in Sections [5] and [6j 
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4. HAAR RANDOM ISOMETRIES VERSUS BLOCK RANDOM ISOMETRIES 

In this section we perform a very basic analysis of the two probability measures we considered 
earlier on the set of isometries from C n to C k (g) C n . Let us first recall their definitions. The 
easiest way to introduce these measures is via an image measure construction, starting from the 
Haar measure on unitary groups. Recall that both these objects are probability measures on 
the set of isometries 

Isom(C n , C kn ) = {V G M knxn (C) | V*V = l n }. 

The Haar measure /iHaar is the easiest to define: it is the image measure of the Haar measure 
on the unitary group lA{kn) via the truncation operation which erases the last {k — l)n columns 
of a fen x kn matrix. In other words, if P G Mfc nX n(C) is the truncation operator Pij = 5ij and 
U G U(kn) is a Haar-distributed unitary matrix, then UP ~ //Haar- 

An isometry V G Isom(C n , C kn ) has distribution fi w if it is obtained by stacking k independent 
Haar unitary matrices one on top of the other, with weights ^/wi: 

k 

(19) V = '£ t y/Wie i ®U®. 

i=i 

The probability distribution [i w can also be seen as the image measure of the product of k Haar 
measures on U(n) via the weighted stacking procedure described above. 

We gather next some basic properties of these measures, whose proofs are left to the reader. 

Proposition 4.1. The measures \xnaar end Ubiock have the following invariance properties: 

(1) IfV ~ HHaar o,nd U G U(kn), U 1 G U(n) are fixed unitary matrices, then UVU' ~ fiHaar- 

(2) IfV~[i w and U' G U(n) is a fixed unitary matrix, then VU' ~ fi w . 

(3) IfV^[i w and U G K(kn) is such that 

k 

i=i 

where a G is a permutation that leaves the vector w invariant and U^> G U(n) are 
unitary matrices, then UV ~ [a w . 

One can easily discriminate statistically between the two measures by computing moments 
or covariances for different matrix entries. In what follows, for the sake of simplicity, we shall 
consider the equiprobability vector w* = (1/k, . . . , l/k). Let us start by computing the moments 
of a single matrix element. It is well known (see |28| ) that, for a Haar unitary matrix U of size 
n, one has 



It follows that, for the two ensembles we consider, we have 

EHaarl^llI 2 ^ 



kn + p — 1 
kn — 1 



^j^ii*=*-'( n ^7 1 )' 



Although the above expressions agree at p = 1, they are different at p = 2, showing a 
statistical difference between the two ensembles. 

More striking examples come from covariance computations: matrix elements Vu and Vi in +i 
are independent under fj, WiK , while this is obviously not true for //Haar, see [281 Proposition 4.2.3]. 
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5. Product of conjugate channels with bounded output dimension 



We start by representing random unitary channels (and the corresponding complementary 
channel) in the graphical formalism we introduced. In terms of the random block-Haar map V, 
the channel has the following form: 



*P0 



i=i 



WiUiXU* = Tv k ((VW ® l n )VXV*(VW ® I n )) 



where X E M n (C) is the input matrix and V = (VW <8>I n )V is the isometry in the Steinespring 
picture. As before, we define the weighting matrix W = diag(wi, W2, ■ ■ ■ Wk) and V is obtained 
by stacking the unitary matrices Ui one on top of the other, as in equation (|19p . Graphically, the 



diagram corresponding to the channel $ is presented in Figure 4(a) , whereas the complementary 
channel is depicted in Figure |4(b)| 

Next, we want to describe the limiting output eigenvalues of a fixed input state going through 
a random unitary channel. We are interested in the following k 2 x k 2 random matrix: 



(20) 

Here, ip n E C n (g) C n is a fixed input vector for each n; notice that we are considering only 
rank-one inputs, since these states are known to yield minimal entropy outputs. 



c $(X) 



$ C (X) 



w 



V 

d 'm — d x 



(a) 



V 



X 



V" 



b 



(b) 



Figure 4. A random unitary channel and its complementary. 

To represent the input ip n in the graphical calculus, we add A n and ^4* boxes on the wires of 
the Bell input, as in Figure [5j 



r 



Figure 5. Generalized Bell states are used as inputs. 
Algebraically, we consider a sequence of inputs 

n 

(21) ip n = ^2 aijei (g> ej, dij € C. 

In the matrix version, this reads 



(22) 



A n — ^ ^ ® &j- 

»J=1 



We require the normalization relation ||^ n || = 1, which is equivalent to 
(23) Tr[A n A* n ] = 1. 
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The input vectors ip n generalize Bell states 

1 - 

which correspond to the trivial choice A n = I n /y/n. 

5.1. Calculation of the limiting eigenvalues. In order to define well-behaved inputs, as in 
|12j . we introduce two assumptions on the asymptotic behavior of the sequence of input states 
-An- 

Assumption 1: 



for some m G C. Note that a similar relation holds for A* n , with m replacing m. Note that one 
has 

m = lim (i/j n ,(p n ), 



n— >oo 
|2\ 



so that one can say that m (or rather \m\ ) measures the overlap between the input state ip n 
and the Bell state ip n . 
Assumption 2: 

(25) PJoo = O ' 



n 

Recall that the empirical eigenvalue distribution of a self-adjoint matrix Z G M fc 2(C) is the 
probability measure 

i=l 

where Ai, . . . , A fc 2 are the eigenvalues of Z. 

Before we state our result, we introduce one more notation, essential to what follows. 

Definition 5.1. Let S : R 2 x R k — > R fc be the function defined by 

S{x,y;w) = spec 1 (Hx(x,y;w)), 

where spec^ denotes the ordered spectrum of a self-adjoint matrix and 

{x + y)wf ifi=j 



(26) Vi,j€[fc], Hv(i,j) 



ywiWj if i / j. 



Theorem 5.2. Under the assumptions (|24p and (|25[) above, the empirical eigenvalue distribu- 
tion of the matrix Z n converges almost surely, as n — > oo, to t/ie probability measure: 



< 27 » F 



i,j=l i=l 



where s = S(l — \m\ , |m| ; w). 

Before we prove this theorem, let us state some of its corollaries and analyze the limit entropy 
of the matrix Z n (which is the entropy of the probability vector appearing in the conclusion of 
the theorem) as a function of the parameters \m\ 2 and Wj. 

We analyze first the "extremal" cases for the weight vector w. 

Corollary 5.3. In the case where the weighting vector is uniform, W = I/k, the S function 
can be evaluated to give 

, . Iml 2 1 — |m| 2 1 — |mp 

(28) Sl = i-L + —^. s l = —^ (2 

This implies that the output state has asymptotically the following eigenvalues: 



12 
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• ^ — h 1 ; with multiplicity one; 

• 1 y;^/^ multiplicity k — l; 

• p-, wzi/i multiplicity k 2 — k. 
The entropy of the probability vector 

( 



I 12 i i 1 2 

m 1 — m 



+ 



1 



1 



k 2 



k 2 



k 2 



V 



Iml 1 1 

^ ; ¥ T 2 

k-1 times k 2 —k times / 

is a decreasing function of \m\ 2 , the asymptotic overlap between the input vector ip n and the 
Bell state ip n . 

Corollary 5.4. In the case where the weighting vector is trivial, w = (1, 0, . . . , 0) the channel 
& is a unitary conjugation and the output matrix Z n is a pure state of null entropy. 

We now turn to the "extreme" values of the parameter m, \m\ = 1 (the input state is, up to 
a phase, a Bell state) and m = (the input state is orthogonal to the Bell state). 

Corollary 5.5. In the case where the input state is equal, up to a phase, to the Bell state, i.e. 
J — 1, the matrix is, up to a constant, a rank one projector and thus 

k 



\m\ 



si 



E 

i=l 



WJ 



S2 



Sk = 0. 



Corollary 5.6. In the case where the input state is orthogonal to the Bell state, i.e. m = 0, 
the matrix Hy, is diagonal and thus 

Si = w 2 , Vi G [k]. 

The limiting eigenvalue vector of the output state Z n is w <£> w and its entropy is thus 

lim H(Z n ) = H{w ®w) = 2H(w). 

n— >oo 

Proof of the Theorem \5.SX The proof uses the moment method and consists of two steps. First, 
we compute the asymptotic moments of the output density matrix Z n and then, by a Borel- 
Cantelli argument, we deduce the almost sure convergence of the spectral distribution and of 
the eigenvalues. 

Step 1: We calculate the limit moments of Z„, using the graphical calculus, see Figure [6l 
Here, o and • correspond to the n-dimensional space, and □ to the /c-dimensional output space. 



V 



A, 



c VW 



v 



At 



v V 



Figure 6. The diagram for the output state Z n . 

In order to compute the p-th moment of the matrix Z n , we use the graphical calculus on 
a diagram obtained by connecting p copies of the diagram in Figure E] in a tracial manner. 
For fixed p G N the Weingarten sum in Theorem 13.111 are indexed by pairs of permutations 
(a,/3) G We label the V and V boxes in the following manner: 1 T ,2 T , ...,p T for the V 
boxes of the first channel (T as "top") and 1 B ,2 B , . . . ,p B for the V boxes of the second channel 
(B as "bottom"). We shall also order the labels as {1 T , 2 T , . . . ,p T , 1 B , 2 B , . . . ,p B } ~ {1, . . . , 2p}. 
A removal r = (a, /3) G 5| of the V and V boxes connects the decorations in the following way: 
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(1) the round white decoration o of the i-th l/-block is paired with the round white deco- 
ration of the a(i)-th V block by a wire; 

(2) the round black decoration • of the i-th. U-block is paired with the round black decora- 
tion of the /3(z)-th V block by a wire; 

(3) the square white decoration □ of the i-th U-block is paired with both the square white 
decorations of the a(i)-th and /3(i)-th V blocks by wires. This double pairing is suggested 
by the duplication symbol associated to the square label. 

We also introduce two fixed permutations 7, 5 G e>2 P which represent wires appearing in the 
diagram before the graph expansion. The permutation 7 represents the initial wiring of the □ 
decorations (corresponding to the trace operation) and 6 accounts for the wires between the • 
decorations connecting boxes A or A*. More precisely, for all i, 

(29) 7 (i T ) = (i - if, 7 (i B ) = (i + l) B , and 5(f) = i B , S(i B ) = i T . 

After the removal procedure, for each pair of permutations (a,/3), we obtain a diagram T> a ^ 
consisting of: 

(1) o-loops; n#( a ) 

(2) u-nets; fw(<*,P) 

(3) •-necklaces; /a(/3) 

First, one can easily see that the number of o-loops is exactly ri& a . Next, for /w(a,/3), since 
square labels are connected by wires with the box T's and T*'s, the graph they yield is not a 
collection of loops, but can be more general, where the boxes Ws are "caught in nets" which 
are made of T's and T*'s. The general formula for fw can be found in Lemma 15.71 Finally, 
the contribution of •-necklaces depends on the moments of the matrices A n and is encoded in 
a function /a(/3) (see |12| for a more detailed treatment of a similar situation): 

(30) f A (J3) = H Tr[A 8 ^---A'^M] 

ceCycleCS- 1 ^) 

Here, |c| is the number of elements in c and s C) i . . . , s c u are defined such that 

, , [ 1 if the ith element in the cycle c belongs to T 

(31) s ci = < 

I * if the ith element in the cycle c belongs to B 

Note that the above function /(/3) is well-defined in spite of the ambiguity of s Cj j, because of 
the circular property of the trace. 

Therefore, the Weingarten formula in Theorem 13.111 reads 

(32) ETr[Z£]= £ n# a f w (a, /3)/ A (/3)Wg(n, a, (3). 

a,0£S 2p 

Using the moment growth assumptions (|24p . (|25p for the matrices A n , we get that, for all cycle 
c of (3, we have 

(33) I Tr [A 110 ' 1 ■ ■ ■ A S °.W] | < \\A S ^ ■ ■ ■ || A^l-l- 1 • ||^4 s -.l-i ||i 

|c||-l 

• yfn = n 1- '^/ 2 as n — > 00, 

where the notation f(x) < g(x) means that there exists some constant C > such that f(x) ^ 
Cg(x) for x large enough. The above inequality is the only place where Assumption 2 (see (|25j)) 
is used. Hence it yields the following asymptotic bound for the factor /(/3); 

(34) \f(P)\<n#^-P. 

Using the equivalent for the (modified) Weingarten function in Proposition 13.9^ we get (note 
that the factors depending on the fixed parameter k are hidden in the < notation) 

(35) ETt[ZP]< ^2 n* a n*^ l5 '>- p n- 2p -^ 1 ^ as n -> 00. 

a,/3eS 2p 
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The power of n in the RHS of (|35p is bounded by using Lemma 13.51 as follows. 

(36) 2p- \a\ +p- |/3 _1 5| - 2p- \a~ l f3\ = p-(\a\ + ja" 1 /?! + IP^^l) < 
Here, equality holds if and only if id — > a — > (3 — > 5 is a geodesic: 

(37) a = J] n, p = n n 

where t% = (i T , i B ) and A C B C {1, . . . , p}; we refer to [14] for a proof of this fact. Importantly, 
for such geodesic permutations /?, the following asymptotic behavior follows form the first 
assumption on the growth of the trace of A n ~. 

(38) f A (f3) = {n\m\ 2 ) m +0{n- 2 ). 

Note that \/3\ = \B\. This implies that the power of n in (|32p in fact becomes for all the a,f3 
which satisfy the geodesic condition id — > a — )■ /3 — > 5: 

(39) #a + \B\ -2p- = 2p - \A\ + \B\ - 2p - \B \ A\ = 

Here, \B\ = \B\A\ + \A\. 

Hence, by using (jlOp . (jlip and (j38j) . we have the following approximation on ETr[Z^] (note 
that the estimate on the error order is not necessary here but will be so in Step 2): 

(40) E[Z£] = (l + 0(n- 2 )) £ f w (a^)\m\ 2 ^ ( _i)l^\A| 

id— s-a— 5-/3 — 

For the above error term, note that Lemma [3.51 implies that permutations (a,/3) G x S2 P off 
the geodesic make the power of n less by two or more; only even powers are allowed. 
Using Lemma 15.81 we can further process the moment expression 

(41) lim ETi[ZP] = {Tr[W p }) 2 

n— ¥oo 

\A\ 

+ Y, Tr 1 W2P ] IH 2|B| (-1) |B| + E Il Tr [W 2 ^-^] |m| 2 l B l(-l)^l. 
A=0,B^0 0^ACB i=l 



(*) (<?) 
The multinomial identities: 

^ x^ = {l + x) p 
0CAC{l,...,p} 
^ x |A| y |B\A| = (1 + x + y)P _ 

0CACBC{l,...,p} 

give further calculations 

(42) (*) = Tr[^ 2p ] [(1 - |m| 2 ) P - l] . 

Set C = B \ A such that C C C [p\ \ A and use Lemma 15.91 with x = 1 — \m\ 2 and y = \m\ 
to get 

(43) (9) = [W 2 ^-^] ]T | m | 2 (|C|+l^l)(_i)|C| 

A^0«=1 0CCC[p]\A 

(44) = ^n Tr [^ 2(ai+i ~ a!) ] |m|2|A| ( i " |m|2 )^ |j4 ' 

A^0 i=l 
k 

(45) = s? i - ( X - l m ! 2 ) P Tr ' 

i=l 
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where, by Lemma 15.91 and Definition l5.il the vector s £ M. k is given by s = S(l — \m\ 2 , \m\ 2 ;w). 
Therefore, whole, 

(46) lim ETr[Z p ] = (Tr[W p }) 2 + (*) + (9) 

(47) = (Tr[W p }) 2 -Tr[W 2p ] + Y^ * 



k 

P 

i 

i=l 



(48) = £(^) p + 5>?, 

i,j=l 1=1 

expression in which one can recognize the limiting eigenvalues announced in the theorem. 

Step 2: We now move on to prove the almost sure convergence. Since this part of proof 
is very similar to that of Theorem 6.3 in [14] or Theorem 3.1 in |12| . we only sketch here the 
main ingredients. Using Borel-Cantelli Lemma, it is enough to prove that the covariance series 
converges: 

OO oo 

(49) 5>[(Tr[Z£]-ETr[Z£]) 2 ] = J> [(Tr [Z£]) 2 ] - (ETr [Z£]) 2 < oo 

n=l n=l 

which will imply that for all p ^ 1 

(50) Tr [Z p ] -> flU a.e. as n -> oo. 

Also, note that by Carleman's condition, equation (|46p uniquely determines the measure as 
in dSTJ). 

First, ()40p implies that 

id— s-a— yf)— >8 

On the other hand, we use Theorem EH] to calculate E[(Tr [Z p ])]. In the dia gram we have 
two identical copies of Tr Z%, which amounts to a total of 4p pairs of V and V boxes. As a 
result, removals (a,j3) are defined for a, [3 G S^p. However, importantly those two copies are 
initially separated. Namely, initial wires 7, 5 G S2p ® S^p = S^ p are written as direct sums: 

(52) 7 = 71 ® 72 and 5 = <5i ® 5 2 , 

where the indices 1 and 2 refer to the first or the second group of 2p boxes appearing in the 
diagram. Then, as before, we calculate the power of n, which is 

(53) 2p- (\a\ + \a~ l p\ + \fi~ l 8\) < 0. 

Equality holds if and only if id — > a — > (3 — ¥ 8 is a geodesic. Moreover, this geodesic condition 
implies that a and f3 can be written as 

(54) a = a x ® a 2 and f3 = $\ /3 2 



(51) (ETr[Z p }) 2 = \k- 2p fw(a,P)\m 




Here, pairs (ai,/3\) and (0:2, ^2) are defined as in ([3] 

Therefore, in the diagram all removals which matter as n — > 00 keep those two copies sepa- 
rated. Also, these removals have the following properties: 

fw(a,P) = fw{a*i,Px) x fw(ai 2 ,P2), 

\P\ = \Pi\ + l/3 2 |, \a~ l fi\ = |arVi| + K'&l 
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The first statement says that fw can be calculated for the each copy independently. Then, for 
the same reasons as before, we get an approximation with the error of order 1/n 2 : 

(55) E(TvlZ?]) 2 = 0(n- 2 ) 

+ k- 4 ? J]) [Av(ai,A)Ay(a a ,A)|m| a ^l +a "AI(-l)K 1 AI+K- 1 Al' 

id— ><5 
id— >a 2 — yffa— >$ 

Finally, we see from (|5ip and (|55p that 

(56) E [(Tr [Z%]) 2 ] - (ETr [Z p ]) 2 = 0{n- 2 ) 

which proves (|49p . and finalizes the proof of the theorem. □ 

The following rather technical lemmas are needed in the proof of the result above when 
dealing with nets containing W boxes. 

Lemma 5.7. For fixed permutations a, f3 6 S2 P 

(57) fw(a,P)= H Triyl"! 

feG(7 _1 aV7- 1 /3) 

Here, we understand the notation j~ 1 a v'7~ 1 /3 in terms of partitions: both permutations 7 _1 a 
and 7 _1 /3 naturally induce partitions on {1, . . . , 2p}, and V is the join operation on the poset of 
(possibly crossing) partitions. Also, b and \b\ stand for a block of a partition and its cardinality, 
respectively. 



Proof. First, we observe that because of the cyclic structure of Z n we can put two yPF-boxes 
into one VF-box and associate it to the neighboring [/-box. Next, we claim that since W is a 
diagonal matrix, we can slide W along wires within the net. Indeed, we can show it algebraically: 

(58) J2(e* W) (g) a fg) a = Wie* <g> a ® a = ^ 4 ® (We») ® e; = ^ e i ® e * ® (Wei). 

i i i i 

Finally, by using the above fact, we collect the W matrices together. However, the rest which is 
composed of T's and T*'s can be contracted to a point. Therefore, each connected component 
b in 7 _1 a V 7~ 1 /3 gives the factor TtW 1 6 I. □ 

The above general formula must be studied in details in order to complete the proof of 
Theorem l5.2[ We need precise values of fwi®-, P)) when a, f3 lie on the geodesic id — > a — > — > S. 

Lemma 5.8. Suppose permutations a, (3 lie on the geodesic id — > a — > /3 — > 5. As in ([57|) . they 
admit decompositions as products of disjoint transpositions indexed by subsets A C B. Then: 

• When A = we have 

'[Tr(W p )} 2 ifB = <H 
Tr (W 2 'P) ifB^Hi. 

• When A = {a% < ■ ■ ■ < a^} ^ ^ we have 

lA{ r 

(60) f w (a,/3) = l[Tv[w 2 ^-^ . 

i=i 

Here, we understand Omi+i — ol\a\ t° be equal to p + a\ — ol\m. 

Proof. First, we consider the case A = 0. If B = 0, the net associated to W is composed of two 
cycles containing each p W boxes (the top and the bottom cycles). However, if B is nonempty, 
the top and bottom cycles become connected and one obtains a large cycle of length 2p. In 
conclusion, these cases yield respectively /w(id, id) = [Tr(VF p )] 2 and fw(id,/3) = Tx(W 2p ), by 
using Lemma 15.71 

Next, let us assume that A = {ai, . . . , aui} ^ 0. The structure of a and [3 as in ([37|) with 
A C B implies that the connected components of the net are determined by indices i £ [p] such 



(59) fw(}d,p) 
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that i G A and i £ B, see Figure [7] for a proof. Since the first condition implies the second 
and each such connected components carries 2(a^ + i — a^j boxes W, one obtains the announced 
formula. 



-q<= 



V 



V 



z>b— c w 1 



i A, i B 
(a) 



— L ^ 



V i 

~z>b- 



-q<— 



V 



V 



i £ A,i e B 
(b) 



V 



-qc±J 

! v 



V i 



V 



i e A, i e B 
(c) 



Figure 7. The connected components of the net containing boxes corre- 
sponding to a pair of geodesic permutations a, (5 is determined by indices i £ A. 
In the picture, only the last case induces a "cut" in the cycles, creating additional 
connected components. 



□ 

The following result is needed to simplify the formulas in the lemma above. It can be, 
however, interesting on its own, from a combinatorial perspective. 

Lemma 5.9. For a real diagonal matrix W = diag(u/i, . . . , w^) and real numbers x,y £ R, we 
have 

\A\ k 

(61) J2 xP ~ lAl y lAl II Tr [W 2 ^- a ^\ =Y,^-x p Tr [W 2p ] , 

0^AC[p] i=l i=l 

where A = {ai < ••• < a|A|}> ( a |A|+i — a \A\) should be understood as (p + a\ — cl\a\) and 
s = S(x,y;w) is as in Definition \5.1[ 

Proof. To the diagonal operator W = Y^i=i w i e i e i> we associate the vector 

k 

R k g> R k 3 W = Wiet <g> a 

i=l 
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and Pyy £ Mfc2(C), the orthogonal projection on W. The idea of the proof is to consider the 
self-adjoint operator 

H = xW ® W + yP w 

and to expand Ti(H p ). Since xW <g> W and yPyy do not commute in general, one has to 
consider general words in these two matrices. Such words can be indexed by the positions 
A = { ai < ■■■ < where P^ appears in the word; let us call Wa the word corresponding 
to a subset AC[p], Oviously, one has Tr(W ) = x p [Ti{W p )] 2 , and thus 

Tr{H p ) = x p [Tv(W p )} 2 + ^ Tt(Wa). 

Using the graphcal notation, it follows from Figure [51 that the general term Tr(WA) factorizes 
along the intervals defined by the set A and one has, for all A ^ 0, 

lAl r 

(62) Tr(W A ) = xP-WyW JjTr I w^ * 1 " *) . 

i=l 



w 



(a) 



(b) 




Figure 8. Diagrams for W (8> W, P^, and for an interval inside a word Wa- 



More precisely, each conected component in the diagram for Wa corresponds to an interval 
+i— ai and, for such a trace (see Figure 8(c) ), one has a contribution of x ai+1 ~ ai ~ 1 yTi [w 2 ( ai+1 ~ ai ^~\ 



Multiplying all these contribution gives ([62]) . 
It follows that 



(63) 



0^AC[p] 



p-I^LI^I 



\A\ 



II*[ 



Tr(tf p ) -x p [Tr(W p )f 



One can further simplify this equation by using the explicit form of the operator H. Indeed, 
either by using the graphical notation or by simple algebra, the action of H on some basis 
vectors can be easily computed as follows. For i ^ j, it is obvious that P^e, g) ej = and thus 
Hei <X> ej = xwiWjei <X> ej, proving that ej ® ej are eigenvectors of H for the eigenvalues xwiWj. 

Let £ = ©^ =1 Mej <8> be the subspace orthogonal to the space spaned by e« ® ej with j 7^ j. 
The restriction of H to S is exactly the operator Hy, defined in (|26p . Definition 15.11 Thus, the 
eigenvalue vector of H% is s = S(x,y;w). We now have computed all the k 2 eigenvalues of H 
and we have 



(64) 



Tr(H p ) = x p {w lWj ) p + 4 



i=l 
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The conlcusion of the lemma follows now easily from (J63J), (|64p . and the following equality 

k 

[Tr(W)] 2 = Y,( w ^r+zZ w ' p - 



i+3 



1=1 



□ 



6. Product of conjugate channels with unbounded output dimension 

In this section, we consider the case where the output dimension grows with the input di- 
mension of the channel, in a linear manner: 

k/n — > c, 

where c is a positive constant that we consider as a parameter of the model. Since both n and 
k grow to infinity, there is no incentive to consider complementary channels, so we focus on the 
output of the original channels 



(65) 



n*](w;), 



where ip n is the generalised Bell state introduced in the previous section. We shall make the 
same assumptions (|24| ) . (|25|) on the growth of the matrices A n appearing in the definition of ip n . 
Moreover, since the number of unitary matrices in <I> grows with n, we introduce the following 
assumptions on the growth of the weight matrices W n . 
Assumption 3: 



(66) 



> 1, lim \ Tr \(kW n ) p ] = t p = [ x p du(x). 

n^oo k J 



where t p are the moments of some given compactly supported measure v. The probability 
measure u, or, equivalently, the moment sequence (t p ) p ^i, are parameters of the model and 
they are fixed. The trace-preserving condition for the channel <E>, Tr(W n ) = 1, implies that 

h = l. 

We first compute the moments of the n 2 x n 2 output matrix Z n . 

Theorem 6.1. Under the assumptions (E 
following asymptotic moments: 

ETr[(cnZ„) 2 ] = 



25]) . and (|66p . the output matrix Z n has the 



(67) 
(68) 



ETr[(cnZ n ) 



t 2 2 + c 2 + tl\m\ A + 0{n~ l ); 
t p 2 \m\ 2p + 0(n~ 1 ), Vp^3. 



Proof. We start by applying the graphical expansion procedure described in Theorem 13.111 to 
the diagram for ETr(Z^), obtained by connecting p copies of the diagram for Z n , displayed in 
Figure O in a tracial way. We obtain a formula which is very close to (|32|) , the only differences 
coming from the fact that we are not using complementary channels. The notation is the same 
as the one in the proof of Theorem [ 



< 








V ' 




r 


V ^ 











[at in 



v 



Figure 9. The diagram for the output state Z, 
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(69) ETr[(cnZ n f] = (any £ n*^^ f w (a, (3)f A (f3)Wg(n,a, (3). 

Let us first upper bound the factor /w(a, /?). The boxes W n appearing in a diagram T> a ^ are 
connected in a net whose connected components are given by the blocks of the partition a V (3. 
In the spirit of Lemma |5, 8\ each such connected component b contributes a factor of Tv(Wn^), 
which, by ([66]) . is equivalent to k^^t^. Hence, 

(70) f w (a,(3)~k#^)-2p TJ i|6| . 

b€aV0 

Replacing k ~ ere and using the bound (|34[) for f A ((3), we can express everything in terms of 
c and n. The power of n appearing in the general term of the moment formula can then be 
bounded by 

P + #(7 _1 «) + #(/?~ ^) - P + #(a V 0) - 2p - 2p - la" 1 /?) 
^2p-(| 7 - 1 a| + |a- 1 /3| + l^ 1 ^l + 
^2p-(| 7 - 1 /3| + ir 1 <5| + l/3|)^0. 

Here, the first inequality holds because 

(71) #(«V/3)^# 
the second one is true by the triangle inequality (Lemma [ 

(72) \ 1 - 1 a\ + \a~ l p\>\ 1 - 1 p\. 

Finally, the last inequality follows from the proof of Theorem 6.8 in [16J, which we recall as a 
lemma. 

Lemma 6.2. For any permutation (3 G S2 P , one has 

(73) | 7 - x /3| + ir^l + |0| > 2p 
equality holding if and only if 

(74) /3 = 

Next, we analyze when the power of n becomes 0. When f3 = 5, the condition (|7ip enforces 
a to be on the geodesic: id — >• a — > (3 = 5. However, with the equality condition for ([72"]) . we 
conclude that a = 5. It is easy to see that f3 = id implies a = id via ([71]) and that /3 = 7 results 
in q = 7 because of the equality condition for (|72j) . One can easily check that for these values 
of (3, the bound (|34[) is saturated from the calculation below. 

Since we identified the dominating terms in the moment equation (j69[) . it is now easy to 
compute the limits; we just plug the following equivalents into (|69j) . 

a = (3 = 5 f w (5, 8) ~ {cn)- p t p 2 f A {5) ~ re p |rre| 2p 

a = /3 = id / w (id,id) = l A(id) = l 

« = /5 = 7 /w(7, 7) ~ (cnf- 2p tl f A ( 7 ) = 1. 

When p ^ 3, we get 

(75) ETr[(cnZ n ) p ] = ^|m| 2p + 0(n _1 ), 
while when p = 2, we obtain 

(76) ETr[(cnZ n ) 2 ] = t 2 , + c 2 + 4\m\ 4 + C^n" 1 ). 

□ 
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The fact that two different behaviours appear in the limits above depending on the value of p, 
can be explained by the presence of eigenvalues on different scales. We first apply the Hayden- 
Winter trick [27\ 125 j to our weighted random unitary channel setting. Note that the following 
proposition applies to any random unitary channel and any input state, at fixed dimension n. 

Proposition 6.3. Consider the output of a rank one input state through a product of conjugated 
random unitary channels 

z n = [*®i](^ n <), 

where the vector ip n £ C n (g) C n is defined as in (|2ip and the channel $ has weights as in 
(TTJ) . Then, one has the following lower bound for the largest eigenvalue of Z n : 

'TrA ' 2 ' 



n U 



Proof. After expanding the sums, one has 

k 

Zn= W i W A U i ® Uj)l/) n lli*{Ui (8) Ujf. 

i,i=l 

If ip n is the vector corresponding to the Bell state 



n 

e,; e;. 



(77) «.-^E' 
one has 

k 

((p n ,Z n (f n ) ^ ^wfiipn, (U{ ® Ui)lp n 1p*(Ui <g> Ui)(f r . 



i=l 

Using the fact that, for all unitary transformations U, one has (U (8) U)<p n = tym we get 
\i{Z n ) ^ (ip n ,Z n ip n ) ^ |(v9 n ,i/; n )| 2 ^w 2 = —^wf. 

i=l n i=l 

□ 

This Theorem 16.31 gives the following lemma as a corollary: 

Lemma 6.4. Ta/ce input states A n satisfying assumption (|24|) and weights W n satisfying the 
scaling (|66p . Then, for all realizations of the random matrix V , the largest eigenvalue of the 
matrix Z n can lower bounded, asymptotically, as follows: 

lim inf Ai (cnZ n ) ^ t2| m | 2 - 

n— loo 

Proof. First, Theorem 16.31 implies that 

lTV/4 I 2 

(78) cn\x(Z n )>- — cn TrPF 2 

n 



Then, we use ([Ml) and (J66J) with p = 2. □ 

Note that this behaviour is consisted with the moments computed in Theorem 16. 11 for p ^ 3. 
In order to investigate the smaller eigenvalues, we analyse the matrix QnZnQni where Qn — 
I n — E n . Here, E = f n fn ^ s the projection on Bell state ([77]) . Precise statements about the 
spectrum will be made later, using Cauchy's interlacing theorem. Before we state our theorem, 
let us recall the definition of a compound free Poisson distribution. 

Compound free Poisson distributions were introduced in [37] by Speicher and the theory was 
further developed in [28J and [34^ Prop. 12.15]. Traditionally, they are defined via a limit 
theorem that mimics the classical Poisson limit theorem. 
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Definition 6.5. Let X be a positive real number and \i a compactly supported probability mea- 
sure. The limit in distribution, as N — > oo of the probability measure 

1 ]S H /i 

NJ U iV^ 

zs called a compound free Poisson distribution of rate X and jump distribution [i , which is 
denoted by tt\^; it has free cumulants given by 

Kp(n\,ij) = Am p (/i), 

where rn p (n) denotes the p-th moment of the probability distribution /z. 

Note that the usual free Poisson (or Marchenko-Pastur) distributions tt c are special cases of 
the above definition, obtained by letting A = c and fj, = 5%. We also introduce the notation for 
the distribution of the square of a random variable: if X has distribution /i, then // x2 is the 
distribution of the random variable X 2 . It follows that the moments of // x2 are 

m p (fi x2 ) = m p (/x) 2 . 

Before going to the the main result of this section, we make a remark on free cumulants. 
The free cumulants and moments of a random variable satisfy the following relation, called the 
moment-cumulant formula: 

(79) m p = Yi K \ b \ 

aeNC(p) beer 

Here, NC(p) is the non-crossing partition and b £ a is a block of the partition a. For more 
details, please see, for example, [MJ Proposition 1.4]. 

Theorem 6.6. The empirical eigenvalue distribution of the matrix {cn) 2 Q n Z n Q n converges in 
moments to a compound free Poisson distribution with rate c 2 and jump distribution u x2 . 

Proof. This proof is similar with the one of Theorem 6.10 in [16] and uses the method of 
moments. 

(80) ^ETv((cn) 2 Q n Z n Q n )P = c 2p n 2p ~ 2 E Tr Y[{I n — E n )Z n 

p 

(81) = cV- 2 ^(-ljl^^ln-l^^lEItlJ^)^ 

g&Fp i=l 



(❖) 

Here, E n = nE n and g G T v = {h : {1, 2, . . . ,p} — > {I n , E n }}. To calculate (0), we need to set 
the natural correspondence of g in S%p, which is denoted by g: When i £ g _1 (I), 

(82) g((i + lf)=i T and g(i B ) = (i + if 
and when % G g^ 1 (E), 

(83) g((i + l) T ) = {i + l) B and g(i B ) = i T 
With this notation we have 

(84) (0)= /w(a,/?)n#«- lQ )/ A (/3)Wg(a- 1 /3) 

a,/3£S 2p 
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Putting everything together and interchanging the two sums, we obtain 



(85) ^ETr((cn) 2 Q n Z n Q n ) p 



(86) = c 2p n 2p ~ 2 1(Sn)l n-l 9 ^ f w (a,(3)n*^ lci > fAffiWgiaT 1 0) 



(87) =c 2 fn 2 f- 2 £ / iy (a,/3)/ A (/3)Wg(a- 1 /3) E^- 1 )' 9 " 1 ^"^" 



(*) 

Importantly, it was shown in [16] that vanishes unless a belongs to the following set: 

£>2p = {vr G : 7r<5 has no fixed point} 
Moreover, it follows also from [16] that for such a £ S2 P one has 

(89) \aS\^p. 

With the estimates ([TO]) and ([31]) for /p^/ and Ja respectively, we have 

^ETT((cn) 2 Q n Z n Q n f < c^n^x 
n z 

J2 ( C n)#( QV «- 2 f^n#(' 3 - 1 ' 5 )-l Q - 1 / 3 l-^ 1 ^)l + #«" la ) (-ljlirWl Mob^-V), 
where t a v/3 is defined multiplicatively over the cycles of the partition a V /3: 

= £|6|- 

&GaV/3 

The power of n in the expression above is bounded by the triangle inequality and (|89p : 

(90) p - 2 + #(q V /3) - - la" 1 /?! - b _1 (^n)| - |^ _1 a| 

(91) ^ 3p- 2 -(|a| + 1/3-^1 + |a- 1 /3| + l5~ 1 (^n)| + |r 1 «|) 

(92) < 3p_2-(|a| + |a- 1 ,5| + | 5 - 1 ( J E n )| + |r 1 a|) 

(93) ^ 2p-2-(|a| + | 5 - 1 ( J B n )| + |r 1 a|) 

(94) ^ 2p-2-( 5 - 1 (^ n ) + |r 1 |) 
In the above each ^ is = respectively if and only if 

(i)#(aV l 8)=#a, (ii) a — >■ /3 — >• 5, (in) |a _1 5| = p (iv) id ->• a -> 5 
We also claim that 

(95) 2p-2-{g-\E n ) + \g- 1 \) ^0 

and ^ becomes = if and only if g = /, called the condition (v). This is from the following fact: 



(96) #g 



2 g = I 

Jg-HE n )\ g^I 

Hence we conclude that the power of n becomes if and only if the conditions (i) to (v) are 
satisfied. 

The condition (v) implies that 5 = 7, which in turn via (iv) forces a to have the following 
structure: 

(97) a = a T ®a B 

Here, id — > a T — > j T and id — > a B — > ^ B where j T = (p T , (p — 1) T , . . . , 1) and 7 s = 
(1 B ,2 B , . . . ,p B ). This result also enforces a similar structure to /3 via (i): 

(98) p = f (3 B 
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Moreover, the special structure of 5 shows that #((a T (Ba B )5) = #(a T ot B ). Here, for the RHS, 
we identify i T = i B . Then, (iii) imposes #(a T a B ) = p, i.e., a T a B = id. 

The condition (iii) with a E S*2p implies that a~ 1 5 is paring. However, note that a moves 
each point within T-group and B-group, and on the other hand, 5 moves each point between 
these groups. So, a _1 <5 has the following structure: 

(99) a^S = Yl n 

ieA 

Here, A is an index set and Tj = (a[, bf). Now, we reinterpret (ii) as 

(100) id -»• a^(3 -> a^5 



which implies, by using (|97|) and (|98|) . that a = (3. 
Hence we now have three conditions: 

(101) (v) g = I, (vi) a = (3, (vii) a = a T © a B with id — > a T — > 7 T and a T a B = id 

and check that these conditions really gives the leading power. The condition (v) implies that 
Ig^ 1 (£n) | = and the conditions (vi) and (vii) lead to 

(102) f A ((3) = f B (a) = 1 

Note that the conditions (v), (vi) and (vii) are in fact neccessary and surficient for the power 
of n to be because /(/?) achieves the bound (jM|) under these conditions. 
Therefore, 

lim \^Tv({cn) 2 Q n Z n Q n y 



n— >oo n 2 



E c2#(QT W- Mob (id) = £ c^ml 

a T in (vii) aeNC(p) 

Y\_c 2 m p (b' x2 ) = m p (-K c 2 M x'z). 



a£NC(p) b£a 



Here, we used the fact that a E S p such that id — > a — > (1,2, ... ,p) corresponds to a non- 
crossing partion of [p\. The last equality above follows from (|79|) and the moment-cumulant 
formula for the compound Poisson distribution of rate c 2 and jump distribution z/ x2 , finalising 
the proof. □ 

Theorem 6.7. The ordered eigenvalues X\ ^ A2 ^ • • • ^ A n 2 of the output matrix Z n have the 
following asymptotic behaviour, as n — )• 00: 

(1) In probability, cn\\ — > t2\m\ 2 . 

(2) The empirical eigenvalue distribution n i l _i Yli=2 °~(cn) 2 \i converges weakly to the com- 
pound free Poisson distribution tt c 2 u x2 of rate c 2 and jump distribution v x2 . 

Here, u is a probability measure defined in (|66|) . 



Corollary 6.8. In the case where the weighting vector is uniform, W n = 1/k, the limiting 
measure is the Dirac mass v = 5\. The largest eigenvalue of Z n behaves as \m\ 2 /(cn) and the 
asymptotic shape of the lower spectrum is given by a free Poisson distribution of parameter c 2 , 
7r c 2 . Then, almost surely, one obtains the following asymptotic behavior for the entropy of the 
output matrix Z n (see [TBI Proposition 6.12] for a proof): 

'21ogn-^+o(l) ifc^l, 

2 



H(Z, ^ - > ~ ~° 2c 



" }21og(cn)-f + o(l) if0<c<l. 



Proof of Theorem \6.7\ First, let Ai ^ . . . ^ X n 2_i be the ordered eigenvalues of Q n Z n Q n . By 
Cauchy's interlacing theorem (for example, see Corollary III. 1.5 of [5]) one has 

(103) A x ^ Ai ^ A 2 ^ . . . ^ A n2 _! ^ A n 2_! ^ A n2 , 

which, together with the conclusion of Theorem 16.61 proves the second statement. 
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Next, we prove the first part. Since c| Tr^4 n | 2 Tr W 2 — > t2\m\ 2 as n — > oo, we define 



(104) 

Then, Proposition 16.31 implies that 
(105) < - mXl 



= \c\TrA n \ 2 TrW 2 - t 2 \m\ 21 



1 < 



(cnXif 



\ 2 -e, 



t 2 \m\ 2 - e n ~"~ (t 2 \rn : 
After taking expectations and using Theorem 16.11 we obtain 

(106) 1 - 1 ! 



Tr[(cnZ n ) 3 } 
(t 2 \m\ 2 - E n f 



t 2 \m\ 2 - e n " (t 2 \m\ 2 - e n ) 3 
Then, there exists a sequence of positive numbers {e' n } such that e' n — > 

(107) E[cnAi] < t 2 \m\ 2 + e' n 
Markov's inequality implies that for any 6 > 0, 

(108) P(cnAi -t 2 \m\ 2 > J) ^ P(cnAi - i 2 H 2 + e n ^ 5) 

/ irin s ^ E[cnAi -t 2 |"i| 2 + e n ] _ + e n 

(109) < - _— ^— 

The other bound is obvious from Lemma [6 



as n — > oo, and 



□ 



7. Conclusions and final remarks 

In this final section, we would like to compare the results obtained in the current paper with 
similar results for non unit-preserving random quantum channels studied in [14} [12] . 

Before going into details, let us first note that the weights appearing in the definition of 
random unitary channels (pQ) represent more parameters that can be chosen to one's convenience. 
This explains why the limiting objects in the current paper are more general than the ones in 

[nunaE]. 

Let us first analyze the case where k is fixed. We shall compare the results in Theorem 15.21 
to the ones in |12[ Theorem 3.1]. To do this, we must first make sure that the channels we 
compare have the same input and output spaces. We must thus enforce the condition t = 1/k 
in [1'2\ Theorem 3.1]. We find that the k 2 limiting eigenvalues are different in the two cases: 



X RC (Z n 



x RUC (z n ) 



\ 



\m\ ^ 1 ///|~ .1 



1 



m 



J*2 i-*3 ' 2-»2 Z-»3 ' ' Z"2 

n n. A. n A. n 



fc 2 — 1 times 



\m\ 2 1 — \m\ 2 1 — \m\ 2 
—r- + 



\ 



1 — \m\ 1 



k 2 



k 2 



V 



fc— 1 times 



1 

k 2 ' k 2 ' ' k 2 

fc 2 — fc times / 



Finally, let us compare output spectra in the regime k/n — > c. We shall make the same 
assumption t = 1/k for the results of |16[ Theorem 6.11] and also we shall consider usual Bell 
states, which imposes \m\ 2 = 1 for Theorem 16.71 In the case where the coefficients w are "flat", 
i.e. v = 1 in Theorem 16. 7\ the results are identical for random channels and random unitary 
channels: the largest eigenvalue behaves like l/(cn) and the lower spectrum has a limiting shape 
7r c 2 on the scale 1/ra 2 . 

In order to understand whether different weights for the unitary operators in ([1]) are more 
interesting for the purpose of finding counterexamples to additivity relations, one needs to 
understand how the minimal output entropy behaves for a single copy of such channels - this is 
the subject of future work [13] . 
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